arXiv: 1505.00407vl [astro-ph.SR] 3 May 2015 


Solar Physics 
DOI: 10.1007/' 


Coronal sources and in situ properties of the solar 
winds sampled by ACE during 1999-2008 

Hui Fu 1 • Bo Li 1 • Xing Li 2 • 

Zhenghua Huang 1 • Chaozhou Mou 1 • 

Fangran Jiao 1 • Lidong Xia 1 * 


© Springer •••• 


Abstract We identify the coronal sources of the solar winds sampled by the ACE 
spacecraft during 1999-2008, and examine the in situ solar wind properties as a 
function of wind sources. The standard two-step mapping technique is adopted 
to establish the photospheric footpoints of the magnetic flux tubes along which 
the ACE winds flow. The footpoints are then placed in the context of EIT 284 A 
images and photospheric magnetograms, allowing us to categorize the sources 
into four groups: coronal holes (CHs), active regions (ARs), the quiet Sun (QS), 
and “Undefined”. This practice also enables us to establish the response to solar 
activity of the fractions occupied by each kind of solar winds, and of their speeds 
and 0 7+ /0 6+ ratios measured in situ. We find that during the maximum phase, 
the majority of ACE winds originate from ARs. During the declining phase, 
CHs and ARs are equally important contributors to the ACE solar winds. The 
QS contribution increases with decreasing solar activity, and maximizes in the 
minimum phase when QS appear to be the primary supplier of the ACE winds. 
With decreasing activity, the winds from all sources tend to become cooler, as 
represented by the increasingly low 0' + /0 6+ ratios. On the other hand, during 
each activity phase, the AR winds tend to be the slowest and associated with 
the highest 0' + /0 6+ ratios, and the CH winds correspond to the other extreme, 
with the QS winds lying in between. Applying the same analysis method to the 
slow winds only, here defined as the winds with speeds lower than 500 km s^ 1 , we 
find basically the same overall behavior, as far as the contributions of individual 
groups of sources are concerned. This statistical study indicates that QS regions 
are an important source of the solar wind during the minimum phase. 
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1. Introduction 


Identifying the source regions of the solar wind is important both as a fundamen¬ 
tal issue in solar physics (Antiochos et at , l2012l) and from the space environment 
perspective (e.g., Luhmann et al., 12002 . and references therein). This practice 
dates back t o the era when the sola r win d was first measure d (Sn yder and 
Neugebauer, 190(1 : Nolte and Roelof, 1973; Neugebauer et al., Il998l . see also 
IPolcttol 120131 for a historic overview). With the solar wind data accumulated 
throughout several solar activity cycles in both near-ecliptic and polar orbits, 
scenarios have emerged as to how the solar wind sources evolve with solar ac¬ 
tivity. This concerns not only the solar winds sampled by individual spacec raft 
but also the solar winds throughout the heliosphere (Luhmann et al.. 120021 ). 

Traditionally, the studies on the solar wind sources start with categorizing the 
winds into the fast (with proton speed s v over, say, 500 kms -1 ) and slow ones 
(v < 500 kms -1 ) (e.g., Schwenn, 2006l l. Regarding the fast solar wind (FSW), 
the coronal sour ce is g enerally accep ted to be coronal ho les (e .g., Krieger, Timo¬ 
thy, and Roelof. fl 97.'ll : Zirker, 19771 : Gosling and Pizzo. f l999 ). Tr acin g the wind 
sampled by Pioneer VI and Vela, Krieger, Timothy, and Roelof ( 197 3d were the 
first to associate the FSW with a coronal hole. Then Zirker (1977) suggested that 
all coronal holes are sources of the FSW. Using the SOHO/SUMER data, the 
outflows at the b ase of polar (Hassler et al., 1999t l and equatorial (Xia, Marsch, 
and Curdt. 120031 1 coronal holes were measured, with the results supporting the 
notion that the FSW originates in coronal funnels (Tu et al, 1200 5!) . On the 
other hand, while examining the ACE and Ulysses data for four Carrington 
rotations during the Cycle 23 maximum, Neugebauer et al. ( 20021 1 concluded 
that a fraction of the FSW originate also from active regions. 

The sources of the slow solar wind (SSW) are substantially more complex. 
While there exists the consensus that the SSWs are associated with coronal 
streamers, debates remain as to exactly where in or around streamers the SSWs 
origin ate. The scenario proposed for solar minimum conditions by Wang et al. 

( 19981 ) suggests that there are two kinds of SSWs, with one originating from 
streamer stalks and the other from just inside coronal holes and immediately 
adjacent to streamers. The former source is con sistent with the outmoving plas- 
moids found by SOHO/LASCO (Sheeley et al.. 1997 : Wang et al. Il998 ). while 
the la tter s ource is corroborated by the SOHO/UVCS measurements (Abbo 
et al., 20 id ) and consistent with the established inverse correl ation of the flow 
tube expansion with the solar wind speed (Wang and Sheeley, 1990l ). However, 
even at solar minimum, this scenario remains to be complemented with the 
expected source of the SSWs from inside streamers, either via direct flow of 
the p lasma from the magnetically open fields in streamer cores (Noci et al., 
19971 ) or via the evaporation of pl asmas from the magnetic arcades in streamer 


helmets ( Suess et all 1 19991 . also in et all . 120051) . Besides, using the method 
of interplanetary scintillation (IPS) tomographic analysis, Kojima et al. ( 19991 ) 
found that yet another SSW source is the unipolar regions in the vicinity of 
active regions (ARs). A further and more direct study associating the SSW with 
ARs comes with Hinode X-ray and EUV spectral observations, where the edge 
of ARs was shown to host persistent up flows with speeds reaching 100 km s -1 
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(Harr a et al .. 120081 ). which may account for up to 1/4 of the in situ SSW (Sakao 
et al., l2007li provided that these upflows e ventually turn into ou tflows. Indeed, 
van DriehGesztel^^^aL j|201^)_(aJso_seeJCulhane_^^_aLl 1 j20lJ_aji^Mandrini 


et al., 20141 1 showed that these upflows may access coronal magnetic fields that 
open into interplanetary space. In addition, using X-ray high temporal-spatial 
resolution images, Subramanian, Madjarska, and Doyle (2 010h found that the 
magnetic reconnection of co-spatial open and closed magnetic field lines at 
coronal hole boundaries creates the necessary conditions for plasmas to flow 
to large distances. This provides an explanation for largely-blue-shifted events 
observed with EIS/Hinode (Madjarska et al ., 2012 /. indicating these plasma 
outflows are also a possible SSW source. Comparing th e rem ote sensing and 
in situ measurements, Feldman, Landi, and Schwadron ( 2005t l suggested that 
the SSW may also arise from the quiet Sun. When it comes to solar maximum 
conditions, SSWs are found to originate from small coronal hole s and active 
regions wher e ope n magnetic field lines exist (Neugeba uer et al., 120021 ; Wang 
and Shee ley, 2003 : Liewer, Neugebauer, and Z urbuchen. [2004 : Ko et al., 2006 : 
Schwenn, 2006 : Wang, Ko, and Grappin. l2009t l 

While the identified coronal sources of the solar wind are diverse, there seem 
to be an agreement on the approaches behind the identification procedure. First, 
unlike the solar wind speed itself, ionic charge states, especially those of oxygen 
and carbon, are suggested to be a telltale signature of the wind sources. Take 
oxygen for example. The abundance ratio 0 7+ /0 6+ measured in the in situ 
solar wind is generally accepted to reflect the electron temperature in the coronal 
sources, given that it does not vary with distance beyond a frac tion o f a solar ra¬ 
dius ab ove t he solar surf ace (O wocki, Holzer , and Hundhausen, 19831 : Biiergi and 
Geiss, 1986 : Hefti et al., 200fll : Landi et al., 2012bl l. Now that the temperatures 
are different in different coronal regions, a comparison of the in situ charge states 
then allows one to ass ociate the in situ wind with a partic ular cor onal source 
(e.g., Zurbuchen et al, 2000l : Zurb uchen , 2001 : Landi et al., 2012al l. With this 
spirit, Zhao, Zurbuchen, and Fisk ( 2009l f divided the non-transient solar winds 
into two categories: those from coronal holes (CH winds) and those from outside 
coronal holes (non-CH winds) with 0 7+ /0 6+ values lower and higher than 0.145, 
respectively. As a result, about 42% of the ecliptic solar wind was found to be 
of non-CH origin during 1998-2008. Second, a model of coronal magnetic field is 
often indispensable. For this purpose, while sophistica ted M agnetohydrodynamic 
(MHD) models are sometimes adopted (Abbo et al.. [2010 1. the potential-field- 
source-surface model (PFSS) and its variants have been in much wider use. On 
the one hand, this practice established the long-term trend of the wind speed 
being in versel y correlated with the lateral expansion of the flow tubes (Wang and 
Sheelev. fl990h . On the other hand, applying the PFSS mode l with an archive of 
the synoptic magnetogram data leads Luhmann et al. (12002 1 to the distribution 
of sources of the heliospheric solar wind as a functio n of so lar activity for nearly 
three activity cycles. In particular, Luhmann et al. 1 20021 1 found that although 
polar coronal holes exist for more than 80% of a solar cycle, they contribute 
to the ecliptic solar winds significantly only during half of a cycle. During the 
other half of a cycle, the near-ecliptic winds originate from mid- and low-latitude 
sources instead. 
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Given the diversity of the wind sources and the complexity of the activity- 
dependence of these sources during a solar activity cycle, the present study is 
intended to examine, in a statistical manner, the fractions taken up by the in 
situ solar winds from various sources from the activity maximum to minimum in 
Cycle 23. To this end, we start with the in situ wind speed measur ements, and 
adopt the standard two-step mappin g proc edure (Neugebauer et al.. 119981 . 2002; 
Liewer, Neugebauer, and Zurbuchen, 20041 1 to trace the winds to their footpoints 
at the solar surface. We then examine the corresponding coronal images recorded 
by SOHO/EIT as well as photospheric magnetograms, and ask the question 
where the footpoints are located: are they located in a coronal hole (CH), an 
active region (AR), or the quiet Sun (QS)? The solar winds are therefore grouped 
accordingly, enabling us to address the question how their in situ properties differ 
and evolve with different activity levels. 

Our study differs from previous studies with similar objectives or similar 
approaches in the following aspects. First, the approach combining a footpoint 
tracing method with the context of c orona l images follows closely the one in 
Liewer, Neu gebauer, an d Zurbuchen (120041) . which is in turn built on Neuge¬ 
bauer et al. ( 1998 . 200^1 where the imaging data were not used. However, while 
Neugebauer et al. (|l998l ) focusedon_^he_C^ule_^2^23jninimum iS _an(^Liewer, 
Neugebauer, and Zurbuchen ( 20041 . also iNeugebauer et alf 20021 ) concerned the 
Cycle 23 maximum, we examine the solar wind dataset that spans the interval 
from 1999 to 2008 in which the declining phase and Cycle 23-24 minimum are 
included. Furthermore, the solar winds from sources other than CHs and ARs 
are paid attention to, and are classified as the QS winds. Second, both this study 
and the one by Zhao, Zurbuchen, and Fisk ( 20091 ) (ZZF09 hereafter) have similar 
objectives in examining the distribution of wind sources in response to solar ac¬ 
tivity. However, the approach for identifying the sources in this study is different 
from the one by ZZF09 where the 0 7+ /0 6+ values are a primary discriminator. 
We note that, given the uncertainties in both approaches, the results of this 
study are meant not to be contrasted with but rather to complement ZZF09, 
with the hope that new light can be shed on the sources of the near-ecliptic solar 
winds. Third, while both using the PFSS model a nd be ing statistical in nature, 
our study differs from the one by Luhmann et al. ( 20021 ) in that we also employ 
the imaging as well as magnetogram data to classify the so urces instead of using 
the locations relative to the equator as in Luhmann et al. i 2()()2 •. Fourth, given 
the considerable interest in and the complexities associated with the sources of 
the slow solar wind, we will analyze the ACE solar winds in general, and examine 
the slow ones in particular. In Section [21 we describe the data and our method 
of analysis. The results are then given in Section [3] Section [4] summarizes the 
present study, ending with some concluding remarks. 


2. Data and analysis 


The two-step mapping pro cedure used in the present study closely 
one in Neugebauer et al. 1 19981 . 20021 ): Liewer, Neugebauer, and 
( 20041 s ). To initiate the procedure, we use daily averages of the solar 


follows the 
Zurbuchen 
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made with the Solar W ind Electron, Proton, and Alpha Monitor (SWEPAM, 
McComaset aj. ^ 1998t l on board the Advanced Composition Explorer (ACE, 


Stone et all lOOsl h Also used are the daily averages of the abundance ratios 
Q 7 +/Q 6 + recorded by the Solar Wind Ion Composition Spectrometer (SWICS, 
a nd t h e mag n etic fie ld measurements with the Magnetic 
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led by 

Gaaa, 


Field Experiment ('MAG. ISmith et al. 1998). Given that wc are interested in the 
non-transient solar winds, one immediate purpose for using the 0 7+ /0 6+ ratios 
is to eliminate from the ACE dataset those intervals occupied by interplanetary 
Coronal Mass Ejections (ICMEs). To d o this , we adopt the same approach as 
in ZZF09 (see also Richardson and Cane. l2004h whereby we discard the data with 
0 7 +/o 6+ ratios exceeding 6.008exp(—0.00578u), in which v is the wind speed 
in km s _1 . A detailed analysis by ZZF09 shows that this criterion adequately 
separates ICMEs from the non-transient ambient winds, being reliable in 83.2% 
of the cases examined therein. The in situ data used in this study span the years 
between 1999 and 2008, hence encompassing nearly half of the Cycle 23. 

The mapping procedure involves two steps. First, the loci of the solar winds 
are found on the source surface, placed at a heliospheric distance of 2.5 Rq as 
implemented by the coronal magnetic held model. This is done via a ballistic 
approach, whereby the longitude correction due to solar rotation is determined 
by the time for a wind parcel to travel from the source surface to the spacecraft. 
Here a constant wind speed is used, and assumed to be the one measured by 
ACE/SWEPAM. The wind parcel is then traced from the source surface to the 
photosphere by following the magnetic held lines computed by using a PFSS 
model, provided in the PFSS package as part of the Solar Software. Instead 
of usi ng the synoptic magnetograms as was done in e.g., Neugebauer et al. 
( 1998l l. this package uses, as the boundary data, the magnetograms measured 
with SOHO/MDI which are updated every 6 hours. It outputs the magnetic held 
vector on a 39 x 384 x 1 92 grid i n spherical c oordinates inside the source surface 
(for details, see lSchriivcr and Dc Rosal . 20031 1. It should be noted that as implied 
by the mapping procedure, the magnetic polarity at the held line footpoint needs 
to be checked against the one measured in situ. Schrijver and De Rosa (120031 1 
found that during 1997-2001, 83% of footpoint polarities matched the interplan¬ 
etary magnetic held (IMF) measurements at the Earth. In this work, we hnd a 
similar behavior: the footpoint polarities are consistent with what is measured 
by ACE/MAG in 81% of the data from 1999 to 2008. To ensure consistency, we 
do not include in our further analysis those dates when the polarities at the two 
ends of the mapping procedure do not match. Table Q] presents, as a function 
of time, the number of daily samples of the non-transient solar wind (second 
column, labeled “All sources”), which is sub-divided into the counts of the solar 
winds from CHs (third column), ARs (fourth), QS (hfth), and Undefined sources 
(sixth). Given in the parentheses are the numbers that correspond to the cases 
where the magnetic polarities match. In total, during 1999-2008, 2124 samples 
are examined in our further analysis, among which 615 (803, 425) samples are 
associated with CHs (ARs, the QS). One can see that a significant mismatch 
takes place in 2007 and 2008. This is possibly understandable given that, close 
to the cycle minimum, the ACE spacecraft was close to the heliospheric current 
sheet. When mapping the winds to the source surface, a small uncertainty may 
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lead to a wrong polarity. For future reference, Table [2] presents the comparison 
between the footpoint polarity and the in situ one for the slow solar winds with 
speeds lower than 500 km s _1 . 

The footpoints are then placed in the context of photospheric magnetograms 
and the EUV images t aken by the Extreme ultraviolet Imaging Telescope (EIT, 


Pcl|dTpudimcn^_et_nl., 1995 ) onboard SO HO ( Solar and Heliospheric Observa¬ 


tory , Domingo. Fleck, and Polandl . 1995l l. While EIT operates at a number of 
passbands (Fe ix/x 171 A, Fe XII 195 A, Fe XV 284 A and He n 304 A), we 
choose the 284 A one because the images recorded in this passband reflect the 
corona at the highest altitude such that coronal holes are more visible. In this 
passband, EIT takes full-Sun images with a pixel size of 2.6" four times a day. 
For consistency, the field line footpoints are compared with images taken at 
around 13:00 UT on the day corrected for the wind travel time. 

The classification scheme is illustrated in Figure [T] where the EIT images 
(the left column) are overplotted with the footpoint locations represented by the 
red crosses. The photospheric magnetograms, on which our scheme also relies, 
are derived from the PFSS model and given in the right column. The scheme is 
detailed as follows. 

A quantitative approach for identifying coronal hole boundaries is imple¬ 
mented, and the winds that have footpoints located within the hence identified 
coronal holes (CHs) are classified as “CH winds” accord ingly. This ap proach , 
which largely follows that in Krista and Gallagher ( 20091 s ) and Ko et al. ( 2014h . 
is illustrated in Figure[2] If a footpoint is located inside or close to an apparently 
dark area, then a rectangular box (the white box in Figurc[2ji) is chosen to enclose 
this part of the dark region and its surrounding area. An intensity histogram 
is constructed, and a multi-peak distribution is then obvious (Figure [2}:). The 
well-defined minimum between the first two peaks then defines the threshold 
for identifying the CH boundary (see the contours in Figure [2ji, also Figure [2 Jd, 
which is the enlarged version of the part inside the box). On the one hand, this 
scheme enables one to objectively define CH boundaries using the EUV images 
in only one passband; on the other hand, it is not influenced by the variation of 
coronal emissions with solar activity. 

A description of some technical details for implementing this scheme seems 
necessary. In practice, we started with asking the question whether there is a 
dark region close to the traced-back footpoint. By “close”, we mean roughly 
“within 100 arcsecs”. If the answer is Yes, we then draw a rectangular box, 
varying in size but typically a few hundred arcsecs across, which encloses both 
a substantial part of the dark region and its surrounding area. The footpoint is 
always within this box. It turns out as long as the box is sufficiently large, its 
size does not significantly influence what one identifies as CH boundaries, for 
the minima in the different histograms pertinent to different box sizes do not 
differ substantially. If the answer is NO, we visually choose the dark area that 
is the closest to the footpoint, and use the same approach to delineate the CH 
boundary (see e.g., FigureQJll). If there is no large apparent EUV CH altogether, 
then we use the threshold found for some obvious CH one or a few days prior 
to this particular day (An example is shown in Figure [T|d 1, where the CHs near 
the two poles are contaminated so significantly that the minimum between the 


SOLA: ver4.tex; 5 May 2015; 0:32; p. 6 





















ACE solar wind sources 


first two peaks in the intensity histogram can hardly be discerned). The box 
is substantially smaller than the disc size, we nonetheless use the threshold to 
delineate CH boundaries throughout the entire solar disc. A space-dependent 
threshold may be more accurate for mapping CH boundaries on the entire disc, 
but our approach suffices given that our purpose is to examine whether the 
footpoint is located inside a CH. Besides, as illustrated by Figure [5^, while a 
single threshold is adopted, the contours outside the box (the dotted lines) also 
outline CHs rather accurately. 

Our association of a footpoint with an Active Region (AR) or the quiet Sun 
(QS) relies on the magnetic morphology of the photospheric regions embodying 
the identified footpoints. The most obvious features on the photosphere are 
magnetic clusters, which are tentatively named “magnetically concentrated area 
(MCA)”. Intuitively speaking, MCAs correspond to strong magnetic fields. To 
make this definition more objective, the absolute value of the radial component 
of the photospheric magnetic field \B r ^\ computed from the PFSS model is 
used. We experimented with different contour levels, |I3 ri ©|,B, used for outlining 
MCAs presented in Figure [T] and the attached movie. In practice, if \B ri Q\g 
is assigned a value 1.5 — 4 times the mean of l-EA,©!, then MCAs become well 
defined. That the MCA morphology is not sensitive to some given \B t ^q\b, as 
long as it is in the mentioned range, suggests that MCAs have sharp boundaries. 
This is understandable considering that MCAs have a strong spatial gradient 
in |-B rj 0 |. The thus-defined MCAs encompass all the active regions numbered 
by NOAA as provided by solarmonitor Q. However, not all MCA patches cor¬ 
respond to a numbered AR. Many of these turn out to correspond to plages 
wi th magnetic fiel d wea ker than concurrent numbered ARs (see section 3.2 


Zharkova et al.l 120051 ). In the three solar activity phases (solar maximum, 


decline, minimum, see Figure ED the contour levels are chosen to be 2.0, 2.5, and 
3.0 times the mean of |Br,©|, respectively. In the MAX and DEC phases the 
threshold is about 10-20 Gs, which is close to the lower bound (15 Gs) adopted 
by Wang, Ko, and Grappin (2009) to identify slow solar winds from ARs. 

With CH boundaries quantitatively defined, when identifying AR and QS 
sources we need only to concern about the regions outside CHs. An AR source 
is defined when a footpoint is located inside an MCA that is a numbered AR 
by NOAA. Likewise, QS sources are defined when a footpoint is located outside 
any MCA. 

With the present grouping scheme, what is unclassified is then named “Un¬ 
defined”, and corresponds to the case where a footpoint is inside some MCA 
that is not numbered by NOAA. These sources may be associated with a decay¬ 
ing/developing AR, but it is also possible that they are distinct from AR sources 
(see, e.g., 06/25-06/27 2005 in the movie where the source is likely a QS one). 
This is why the word “Undefined” is chosen. 

Such a scheme will not overestimate the counts in the respective groups. First, 
the counts of AR and QS sources are not overestimated, since some footpoints 
deemed “Undefined” may in fact be AR and QS sources. Second, the counts of 


1 http://solarmonitor.org 
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CH winds are not overestimated either, for the current definition of CHs excludes 
a fraction of CHs with overlying bright emissions. In any case, the counts in the 
Undefined group account for only a minor fraction of the samples (11.2%, 9.5%, 
14.5%, 10.4%, 9.3%, 18.5%, 17.9%, 19.5%, 18.0%, and 5.3% for the years 1999 
to 2008, respectively). 

Before proceeding, several remarks on our approach seem in order. The first 
remark is on the reliability of the PFSS model, given its apparently oversim¬ 
plification of imposing a spherical source surface, neglecting volumetric electric 
currents between the source surface and the photosphere, and supposing purely 
radially directed field lines outsid e the source surface. Nonetheless, a detailed 
comparison study by Riley et al. (200d) demonstrated that the magnetic field 
configuration computed by the simple PFSS model agrees well with the one 
found in sophisticated MHD computations, provided that both models are driven 
by the same line-of-sight magnetograms. From the practical point of view, the 
magnetically open regions obtained by the PFSS model well match the c oronal 
hole regions ident ified in, say , the He I 10830 synot pic dia grams ('Levine. Il982 : 
Neugebauer et al., 1998 . 2002 ; Schrijver and De Rosa. [2003 1. A good way to make 
sure that the traced-back footpoints are reasonably accurate is to compare the 
current SolarSoft PFSS results with some other calculations. To address this, we 
randomly chose three Carrington Rotations in the MAX, DEC, and MIN phases, 
and compared our derived footpoints with those derived from the PFSS model 
where the magnetogram input is from the Wilcox Solar Observatory (WSO). 
We found that the fraction of the days when the two different sets of footpoints 
belong to the same open field region is 80.7% for CR 1969, 84% for CR 2005, 
and 82.6% for CR 2054. Nevertheless, let us stress that the fraction that the two 
sets do agree is substantial enough that the statistical study we conduct can be 
deemed reliable. 

Another source of uncertainty may come from the mapping procedure, par¬ 
ticularly in view of the simple ballis tic tre atment in volved in the first s t ep. A s 


demonstrated by Nolte and Roelof ( 1973h (also see Neugebauer et al. . 2002 1 


while the solar wind may experience some acceleration beyond the source sur¬ 
face, this effect may be counter-balanced by the near-Sun corotation. Actually a 
further evidence lending us confidence with this mapping procedure is that, when 
inspecting the footpoints on a consecutive basis (please see online animation 
1 attached to Figure [U, one can see an orderly distribution of the locations 
of footpoints. They stay in a particular group for several days before moving 
to another group. In addition, assuming that the uncertaint y with the source 
longitude at the source surface is ±10° (Nolte and Roelof, 11973 1. we select 
two Carrington rotations in each sub-interval (see Figure [3ja,) and examine how 
well our classification scheme works. This is done by tracing the photosphcric 
footpoint from a locus on the source surface 10° eastward or westward of the 
nominal locus, and then examining whether the footpoint is located in a different 
area in the EIT images. We found that at maximum activity, about 30% of the 
footpoints indeed are associated with an area different from what we identified 
using the nominal locus. During the declining and minimum phases, however, 
this mismatch reduces to < 20% and 10% of the cases examined, respectively. 
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3. Results 


3.1. Sources of the ACE wind between 1999 and 2008 


Having categorized the winds, we then address the question of the percentage 
each kind of wind occupies, and how this evolves in response to solar activity. 
Figure [3Ji presents the monthly average of the smoothed sunspot numbers from 
1999 to 2008. Three sub-intervals, labeled MAX, DEC and MIN, are then defined 
according to the level of solar activity. Figure EJ) presents the percentage of the 
CH (shaded blue), AR (red), QS (green) and Undefined (orange) winds in this 
period. These percentages are yearly averaged values, and add up to unity in 
each year. One can see that in the sub-interval MAX, the QS supplies only a 
small fraction of the winds sampled by ACE (~ 10%), the contributions from 
CHs is ~ 15 — 20%, and more than half (56%) of the ACE winds originates 
from ARs. In the declining phase (2002 - 2006), the contributions from CHs and 
ARs both amount to roughly 34%, and the contributions from ARs (QS) tend 
to decrease (increase). As for the sub-interval MIN (2007 - 2008), the fraction 
of the winds from ARs is only marginal (< 17%), while some 31% comes from 
the CHs, and nearly half of the winds originates from the QS. 

Despite the differences in the approaches for identifying the solar wind sources, 
Figure [3] agrees with Figure 1 in ZZF09 in that there exists a tendency for 
coronal sources other than CHs to contribute significantly to the ACE solar 
winds between 1999 and 2008. Overall, FigurejSJr indicates that non-CH sources 
may be more important than CHs in terms of their mass supply to the solar 
wind. In particular, Figure [3jo indicates that the majority of the near-Earth 
solar wind comes from ARs during the Cycle 23 maximum. This behavior agrees 
with ZZF09, and seems to persist to the Cycle 24 maximu m as indicated by the 
very recent study by Brooks, Ugarte-Urra, and Warren (12015 1. Furthermore, 
Figure [3] indicates that CH winds tend to dominate in the year of 2003, which is 
also in line with ZZF09. What is new in Figure[3]is that combining the imaging as 
well as the magnetogram data allows us to further determine, among the sources 
outside CHs, the fractions of the winds from ARs and the QS. Some apparent 
differences from ZZF09 arise as a result. Around the Cycle 23-24 minimum, 
while ZZF09 indicated that CH winds dominate, our Figure [3)3 suggests that 
the contribution from the QS is more important. We stress that this apparent 
discrepancy stems from the differences in the schemes for classifying the solar 
winds. 

More insights can be gained by examining how the in situ properties of the 
solar winds categorized by our scheme depend on solar activity. We sort the wind 
speeds v into 6 bins uniformly spaced between 200 and 800 km s^ 1 , group the 
0 7+ /0 6+ ratios into 6 bins uniformly spaced between 0.0 and 0.6, then present 
in Figure [I] a contour plot in the v-0 7+ /O e+ space the counts of the winds from 
different sources as labeled. The left, middle, and right columns correspond to 
the intervals MAX, DEC, and MIN, respectively. Consider the interval MAX 
(left column) first. One notices that the majority of the winds corresponds 
to an 0'+/0 6+ ratio larger than 0.145, which we recall is the criterion that 
ZZF09 employed to separate CH winds from non-CH ones. However, a more 
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detailed analysis like ours indicates that not all winds that have 0 7+ /0 6+ ratios 
lower than the nominal value of 0.145 are from CHs. Conversely, winds with 
0 7 +/0 6 + exceeding 0.145 are not necessarily non-CH ones. Now consider the 
years 2007 and 2008, labeled MIN. One can see from Figure 0] (right column) 
that the 0 7+ /0 6+ ratios tend to be low, with the majority being lower than 
0.145, meaning that if categorizing the ACE winds by this threshold, one would 
find that nearly all the winds are from CHs. However, combining the footpoint 
tracing approach with the EUV and magnetic field data, we find that the QS 
is the primary contributor to the ACE winds during this period. Furthermore, 
comparing Figures [4jcl with [4jc3, one can see that the QS winds are distinct from 
the CH winds in that they tend to be substantially slower. To select the proper 
subset of the fast solar wind sampled by ACE that comes from CHs, it would 
be almost unmistakable to choose those with speeds higher than 600 km s _1 
and 0' + /0 6+ lower than, say, 0.05. The contamination from the QS winds 
would be at most marginal, and that from the AR winds would be minimal. 
We note in passi ng that this practice has been successfully employed by Zhao 
and Landi (2014). Regarding the declining phase (middle column of Figure 0|), 
one finds that the possibility of distinguishing between CH winds and non-CH 
winds lies in between the extremes of maximum and minimum conditions. This 
is particularly true in the speed dimension. The CH winds tend to be faster than 
the non-CH ones (mainly from ARs in this case), and the difference between the 
two tends to be more obvious than for the MAX phase, but appears significantly 
less obvious than for the MIN phase. 

The 0' + /0 6+ ratios for the CH winds during the MAX phase (Figure [4jil) 
require some explanation. There appears to be a fraction of the CH winds for 
which the 0 7+ /0 6+ values exceed 0.26. If assuming ionization equilibrium, this 
would correspond to a freeze-in temperature exceeding 1.58 MK (Mazzotta et al., 
199 8). This is beyond the currently accepted electron temperatures derived from 


remo te sensing measur ement s for CHs below 1.6 R 0 (Habbal, Esser, and Arndt, 
1993 : Esser and Edgar. l2000l and references therein). This apparent discrepancy 


is not too worrisome given that this fraction of the CH winds tends to origi¬ 
nate from the boundaries bet ween CHs and ARs, while the measurements made 
by Habbal, Esser, and Arndt (11993 9 pertain to t he reg ion well inside a polar CH. 
Furthermore, as proposed by Esser and Edgar ( 2000 ). the electron distribution 
function may rapidly develop a non-Maxwellian character within the first several 
solar radi i that eventually forms what is measured in situ as the halo electrons 
(Marsch, l2006i ). It is worth noting that this non-Maxwellian character is also 
possible to develop in AR and QS winds. 

The differences in the in situ properties of the winds from different sources 
are further examined in Figure [5j where (a) the wind speed and (b) the oxygen 
charge state ratio are plotted as a function of time. Given by the green, red, 
and blue curves are the parameters of the QS, AR, and CH winds, respectively. 
The standard deviations are given by the error bars for the corresponding values, 
which are slightly displaced from one another for display purposes. An immediate 
impression from Figure [5] is that the CH winds tend to be the fastest, while 
the AR and QS winds have almost the equal speeds. And the 0 7+ /0 6+ ratios 
are lowest (largest) for CH (AR) winds, with the QS winds lying in between. 
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However, the considerable overlap in either the speed or 0' + /0 6+ ranges means 
that neither of the two parameters, on its own, seems to suffice to discriminate 
the wind sources. Regarding the activity-dependence of the parameters, one can 
see from Figure [5ji that the wind speeds in the three categories show a similar 
non-monotonic behavior. Take the CH winds for instance. Their speed start from 
relatively low values (~ 500 km s -1 ) around MAX, rise to some 590 km s -1 in 
2003 before decreasing to around 500 km s -1 in 2004-2006, and then gradually 
increase to some 550 km s -1 toward the MIN phase. Moving on to Figure[SjD, one 
can see that the overall tendency of 0' + /0 6+ ratios in response to solar activity 
is opposite to that of the speed, as would be expected given th e well -established 
inverse correl ation o f the two parameters (Wang and Sheeley, 120031 : Wang, Ko, 
and Grappin, 2009l l. Nonetheless, one can see that the 0 7+ /0 6+ ratios from 
different groups of winds differ more significantly than the speeds do: Note the 
marked difference in the 0 7+ /0 6+ values in the CH winds from those in the AR 
winds in the whole period. It is noteworthy that with decreasing activity, the 
0 7 +/0 6+ ratios in all three types of solar winds tend to decrease, which agrees 
with Lepri, Landi, and Zurbuchen ( 2013l l. The 0 7+ /0 6+ values in AR (CH) 
winds are 0.26 (0.16) during MAX, and decrease to 0.10 (0.05) during MIN. 


3.2. Sources of the ACE slow wind between 1999 and 2008 


Given the considerable interest in understanding the origins of the slow solar 
wind (SSW), it is informative to apply the same practice to the slow winds 
alone. In the present study, a SSW is defined to be the wind with speeds lower 
than 500 km s _1 . Figure [6] examines the time evolution of the fractions of the 
SSWs coming from various sources during 1999-2008. Overall, the impression 
in the MAX and DEC phases is similar to what one finds in Figure [3] where 
the solar winds as a whole were considered. This similarity to Figure [3 Jd is not 
surprising given that, as shown in Figure |4j most of the solar winds is on the 
lower side when the speed is concerned. During the MIN phase, the contribution 
from the QS to the slow wind is even more important than that to the overall 
solar wind. This is also understandable in view of Figure QjH, given that the 
solar winds from CHs are largely fast ones. 

Figure [7] presents (a) the wind speeds and (b) the 0 7+ /0 6+ ratios for the slow 
solar wind as a function of time. As far as the wind speeds are concerned, one can 
see that the speed in a given group does not show a systematic variation with 
solar activity. In addition, there is no clear-cut difference in the speeds of the 
winds from different groups. A stronger temporal variation and a more significant 
difference in different groups of winds he in the 0 7+ /0 6+ values (Figure [TJd) . 
Overall, the 0 7+ /0 6+ values for all the winds show a decrease with decreasing 
solar activity, and they are substantially different for different groups. The dif¬ 
ferences in the 0 7+ /0 6+ values in winds from different sources may be a result 
of the intrinsic difference in the respective source properties, the magnetic field 
strength being the most l ikely one. At any rate, this reinforces the notion raised 
by Antiochos et al. (2012)4 , who suggested that the compositional properties and 
temporal variability serve better in differentiating the wind sources than the 
speeds. 
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The S SW pr operties may be compared with previous studies. Wang, Ko, and 
Grappin (20091) suggested that the slow wind during 1998-2007 mainly contains 
two components: one from small holes located in and around ARs with high 
0' + /0 6+ ratios during maximum, the other from the boundaries of large CHs 
with intermediate 0' + /0 6+ values. Our approach suggests that the majority 
of the former component indeed comes from ARs during maximum. However, 
the latter component may actually come from all the three kinds of sources (see 
Figure |4|. 


4. Conclusion 


The main purpose of this work is to examine, in a statistical sense, the sources 
of the solar wind sampled by ACE during 1999-2008 in general, and those of 
the slow solar wind in particular. To this end, we start with the in situ wind 
speed, and find the photospheric footpoints of the wind pa rcels by em ploying the 
standard two-step mapping procedure (Neugebauer et al., 19981 . 2002 1 whe re the 
Potential Field Source Su rface (PFSS) model (Schatten, Wilcox, and Ness. il969 : 
Altschuler and Newkirk, 19691 ) is used. We then associate the footpoints with 
various areas in the EUV images recorded by EIT in its 284 A passband and 
photospheric magnetograms. With this association we classify the ACE winds 
into three groups: coronal hole (CH), active region (AR), and quite Sun (QS) 
winds. Our main results can be summarized as follows. 


i) During Cycle 23 maximum (years 2000 and 2001), ARs are the main con¬ 
tributor to the ACE winds, the contribution of CHs (QS) is ~ 20% (13%). 
The winds in this interval tend to be slow, and the AR winds correspond to 
substantially higher 0 7+ /0 6+ values than the CH winds. During the declining 
phase, the contributions from CHs and ARs both amount to roughly one 
third. Overall, the fraction of the QS winds in this period is 17%, and tends 
to increase with decreasing activity, accounting for 31% of the winds in 2006. 
During the Cycle 23-24 minimum (2007 and 2008), the contribution of CHs 
(ARs) is about 31% (15%), while the QS contribution is ~ 41%. 

ii) Overall, in each phase of solar activity, the winds from CHs tend to be the 
fastest and associated with the lowest 0 7+ /0 6+ ratios. While both lower 
than CH winds, the speeds of AR and QS winds do not show a substantial 
difference. A slightly more pronounced difference between AR and QS winds 
is seen in their 0' + /0 6+ values, with AR winds tending to be slightly hotter. 
As for the dependence on solar activity of the winds from the same sources, 
overall with decreasing activity the winds tend to have lower 0 7+ /0 6+ ratios. 

iii) The fractions occupied by the slow solar winds from different groups show a 
dependence on solar activity similar to the case where solar winds from all 
speed ranges are considered. This can also be said for the activity dependence 
of 0'+/0 6+ values. During the minimum phase, the QS contribution to the 
slow wind is even more important than its overall contribution, amounting to 
- 47%. 
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Our results suggest that the quiet Sun is an important source of the ACE 
solar winds around the cycle 23-24 minimum. A further study dedicated to 
the examination of the properties of the source region and in situ properties of 
this particular QS wind is needed. Regarding the source regions, such properties 
as the magnetic field strength as well as magnetic topology will be of interest. 
Regarding the in situ properties, the abundances of low first-ionization-potential 
(FIP) elements relative to their photospheric values will be in forma tive (Feld¬ 
man, Landi, and Schwadron, 12005 : Wang, Ko, and Grappin, 2000 1. In addi¬ 
tion, it will be worthwhile looking for direct signatures of outflow in the QS 
by examining the Doppler shifts with the emis sion lines measured with either 
SOH O/SUM ER (e.g., Xia, Marsch, and Curdt. 2003 1 or Hinode/EIS (e.g., Fu 
et al.,[2_014; Brooks, Ugarte-Urra, and Warren. 120151 ). 
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Table 1. Number of daily solar wind samples analyzed in each year. 


Year 

All sources 

CH winds 

AR winds 

QS winds 

Undefined 

1999 

237 (188) 

32 (27) 

148 (117) 

28 (23) 

29 (21) 

2000 

261 (221) 

50 (47) 

145 (125) 

36 (28) 

30 (21) 

2001 

272 (220) 

43 (39) 

155 (119) 

38 (30) 

36 (32) 

2002 

289 (259) 

65 (62) 

152 (135) 

38 (35) 

34 (27) 

2003 

277 (259) 

142 (138) 

71 (65) 

34 (32) 

30 (24) 

2004 

242 (211) 

63 (61) 

94 (81) 

34 (30) 

51 (39) 

2005 

250 (201) 

78 (70) 

70 (61) 

56 (34) 

46 (36) 

2006 

262 (200) 

60 (56) 

55 (43) 

94 (62) 

53 (39) 

2007 

257 (178) 

68 (57) 

36 (30) 

114 (59) 

39 (32) 

2008 

259 (187) 

61 (58) 

29 (27) 

157 (92) 

12 (10) 

Sum 

2606 (2124) 

662 (615) 

955 (803) 

629 (425) 

360 (281) 


Table 2 

. Number of daily slow solar wind samples analyzed 

in each year. 

Year 

All sources 

CH winds 

AR winds 

QS winds 

Undefined 

1999 

188 (141) 

20 (15) 

119 (89) 

23 (18) 

26 (19) 

2000 

204 (167) 

26 (23) 

123 (103) 

29 (23) 

26 (18) 

2001 

245 (194) 

37 (34) 

139 (103) 

37 (29) 

32 (28) 

2002 

238 (208) 

40 (37) 

135 (118) 

32 (29) 

31 (24) 

2003 

111 (96) 

35 (31) 

36 (32) 

20 (18) 

20 (15) 

2004 

194 (163) 

36 (34) 

83 (70) 

28 (24) 

47 (35) 

2005 

172 (132) 

34 (29) 

52 (44) 

47 (28) 

39 (31) 

2006 

205 (146) 

38 (34) 

40 (28) 

81 (50) 

46 (34) 

2007 

187 (116) 

38 (28) 

30 (25) 

91 (42) 

28 (21) 

2008 

169 (102) 

25 (22) 

15 (13) 

122 (61) 

7 (6) 

Sum 

1913 (1465) 

329 (287) 

772 (625) 

510 (322) 

302 (231) 
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Figure 1. Illustration of the classification scheme of the ACE solar winds. The footpoints of 
the solar wind flow tubes are given by the red crosses, which are classified as being associated 
with a coronal hole (the first row), an active region (second), the quiet Sun (third) and some 
undefined source (bottom). The left column presents the EIT 284 A images, while the right 
column gives the corresponding magnetic morphology of the photosphere. The green contours 
outline CH (left column) and Magnetically Concentrated Area (MCA) boundaries (right). An 
animation showing the sources during 1999-2008 is available online. 
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Figure 2. The scheme for outlining coronal holes. Panel (a) presents the EIT 284 image on 
2003 June 11, when the traced-back foot point (the red cross) is located close to a low-latitude 
CH. The white box encloses the region for which the intensity histogram is constructed and 
presented in panel (c), where the solid green line represents the minimum between the two 
peaks, given by the two red dotted lines. This minimum is used as the threshold to delineate CH 
boundaries in (a), where the contours inside (outside) the box are given by the solid (dotted) 
lines. Panel (b) is an enlarged version of the part enclosed by the box in (a). 
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Figure 3. Fractions of the ACE winds with different sources as a function of time. Panel 
(a) shows the temporal evolution of the smoothed monthly sunspot number during 1999-2008, 
which is further divided into the maximum (labeled MAX), declining (DEC) and minimum 
(MIN) phases. Panel (b) gives the percentage of the coronal hole (CH, blue), active region 
(AR, red), quiet Sun (QS, green) and undefined (UN, orange) winds. 
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Figure 4. Dependence on solar activity of the distribution of solar winds from different sources 
in the speed-0 7_ *"/0 6- ’" space. The left (middle, right) column corresponds to the maximum 
(declining, minimum) phase, while the first (second, third) row represents the winds from 
coronal holes (active regions, the quiet Sun). Here the counts of solar wind samples in different 
groups are shown as contour plots with the contours equally spaced in each panel. 
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Figure 5. In situ properties of solar winds from different sources as a function of time during 
1999-2008. Here panels (a) and (b) are for the wind speeds and 0 7 ~*~/0 6_t_ ratios, respectively. 
The interval between 2000 and 2008 is further divided into three activity phases: maximum 
(MAX), declining (DEC) and minimum (MIN). The winds from coronal holes (CHs), active 
regions (ARs) and the quiet Sun (QS) are represented by the blue, red, and green curves, 
respectively. As for the error bars, they represent the standard deviations in each year. 
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Figure 6. Similar to Figure[3]but restricted to the slow wind with speeds less than 500km s~ 1 . 
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